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Abstract 

Based on the study of the simple Abelian Higgs model in 1 + 1 dimensions we will present 
a new method to identify and localize extended instantons. The idea is to measure the 
topological charge on regions somewhat larger than the extended instantons so as to average 
out the ultraviolet fluctuations but without losing the detailed topological information when 
going to the full space. The instanton size and probability density can be directly extracted 
from this analysis. Local dislocations, which can be avoided for fine enough lattices, can be 
reinterpreted as modified boundary conditions producing sectors with net topological charge. 

1 Introduction 

Topological non trivial configurations, such as instantons [1] are expected to play an important 
role in such diverse phenomena as net Baryon number generation [2, 3, 4], the U(l) problem 
[2, 5] or even quark confinement [6, 7]. After the first attempts to calculate topological effects 
with semi classical methods [8, 9], where it is very hard to find the relative weight for the various 
configurations or to go beyond the dilute gas approximation, several works have tried to extract 
directly the topological content from lattice simulations [10, 11, 12, 13]. Up to the present, most 
simulations have been done on fairly coarse lattices where the instanton sizes are of the order of 
the spacing, thus making it hard to distinguish between a physical instanton and a lattice artifact 
(dislocation). Furthermore, the topological charge is usually only measured on the whole volume 
where, as a consequence of the periodic boundary conditions, we can only see a net charge if 
we have dislocations (which per definition are points where the charge gets shifted). In several 
works cooling [15, 14] and improved actions [17, 16] are being used in order to smoothen out 
the instantons and filter the physical ones, leaving some uncertainty on the agreement of the 
topological content with the one in the starting model. 

In this work we want to present an alternative method which extracts directly the topological 
information from the lattice configurations. It consists in measuring the topological charge on 
volumes larger than or of the order of the relevant instantons. Let us first assume that one 
can reach a coupling regime where the gauge configurations are smooth enough so that there 
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are no individual sites where the topological density gets even close to ±1/2. Then one has 
safely configurations without dislocations and if we take periodic boundary conditions the total 
geometrical topological charge is exactly zero. This does not mean that there are no instantons 
in the system, but rather that there is a number of instantons and anti-instantons with total 
charge zero. Furthermore, in the thermodynamic limit of very large volumes V, the topological 
charge probability distribution in some finite sub-volume v should reach a fixed limit P V (Q). 

This probability distribution to get some charge Q in a volume v, which will be obtained 
from the lattice simulation, should have a broad envelope for very large volumes, which scales 
like Pav(Q) = ^ I y/oiPv{Q I y/ot) with av, with distributions centered at integer Q's reaching 
the envelope. As we reduce the volume v to sizes smaller than the instanton, we see only a 
distribution peaked at Q = as we have only pieces of instantons in v or fluctuations. In fact, 
for very small volumes one sees only the ultraviolet fluctuations in the topological charge giving 
a narrow Gaussian distribution (a <C 1). Interestingly, these are topological fluctuations in 
regions of stable topological vacuua and therefore tend to cancel out over fairly short distances. 
Due to their topological character the total dispersion in the charge that they produce only 
grows as yfs with s being the free surface. For larger volumes, of the order or somewhat larger 
than the instantons, we start getting relevant contributions to P V (Q) from instantons, where 
the charge adds coherently to ±1 (if they are fully in v), thus broadening clearly the charge 
distribution in contrast to the fluctuations. Note that instantons extremize locally the action 
and therefore should not be considered as fluctuations themselves. At even larger volumes one 
can do a careful analysis using Poisson distributions of almost free multi-(anti)instantons with 
some average density p and instanton size vj, in order to determine these parameters for a given 
simulation. Then we can check the scaling with v as described above. The general idea is 
therefore to look at volume scales v, where 

a D < vi < v < V 

so that we have mainly supressed the fluctuations and can study the ensemble of (anti)instantons, 
while having effectively almost free spatial boundary conditions. We have tested these ideas 
for the abelian Higgs model in 1 ± 1 dimensions where one has a simple system with non- 
trivial topology. We find indeed, for couplings corresponding to fine lattices, clearly identifiable 
extended instantons over several lattice spacings, while the total charge over the whole V stays 
always zero. Assuming an almost free ensemble of (anti)instantons, we have been able to extract 
an instanton average density p and size, also as function of temperature, T = 1/N T . The density 
p which can also be easily extracted from the second moment of P V (Q) for large v, corresponding 
to the topological susceptibility, falls steeply at higher T as the instantons get squeezed in the 
time direction [24]. We have also studied the evolution of the position of the instantons in 
Monte Carlo time and found that they get created locally as instanton - antiinstanton pairs 
and annihilate again after some time or drift away, to get destroyed later when they encounter 
another pairing. 

As for some models it might get prohibitively expensive to go to such fine lattices where 
the topological density does not get close to ±1/2 for some sites (it might even be impossible 
for cases where the instantons stay almost scale invariant after renormalization) we wanted to 
analyze carefully what happens when at some individual plaquettes (hypercubes in higher D ) 
the geometrical definition for the topological density surpasses ±1/2 before taking the "mod" 
operation, which brings it to the interval ( — 1/2,1/2]. At those points where one had to make 
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the "mod" operation, one introduces "defects" to the topological configuration in such a way 
that the total charge is not zero anymore. We will come to the conclusion, for the abelian Higgs 
model, that these defects only make a difference when there is a net number of them as if we 
had new boundary conditions with this number of instantons, otherwise their effect cancels out. 
In cases with net number N of defects in a total volume v the remaining configuration behaves 
just as if we were in a case without singularities in a sub-volume v with exactly N instantons. 
Furthermore, if we count in Monte Carlo time with which recurrence the configuration stays in 
each topological sector we find a probability distribution very similar to the P V (Q) one without 
singularities. 

In our discussion up to this point we have taken the geometrical definition for the topological 
density which always gives (or an integer for dislocations) in a periodic box. This geometrical 
definition gets fairly cumbersome in non-abelian theories in higher D [19, 20] but should still 
give integer results and our whole idea should work there too. We have also investigated the 
field theoretic definition for the topological density which is easier to implement [21] and gives 
the right naive continuum limit but does have only approximately right addition properties and 
therefore does not give exactly integers on the whole volume V. Nonetheless, on volumes v 
somewhat larger than the instantons as we only need for our analysis, and for not too coarse 
lattices, the field theoretic definition gives good agreement with the geometrical one. 

In the next section we will discuss the abelian Higgs model in 1 + 1 D and its topological 
density definitions. In section three we will choose the parameters in the model in order to have 
effectively a fine lattice. In the fourth section we will identify the instantons in our simulations 
and present a simple model that fits the topological charge distribution. In the fifth section we 
will discuss the role of dislocations and the field theoretic definition for the topological charge 
density. Finally we will discuss our results and propose the use of this method in the non-abelian 
case. 

2 The Model and defining a Topological Density 

Starting point of our analysis is the lattice formulation of the abelian Higgs model in 1 + 1 
dimension with the Euclidean action 

S = -PY, Re(U n ) + i(*n*n ~ I)) 2 - 2« £ Re(<S>* n U n ^ n+ ,) . (1) 

n n ri,/i 

Here & n is the complex Higgs field, U ntfl = exp(i# nj|U ) are the usual links and U n = IinU ntll the 
plaquette variables. (Note that in two dimensions the plaquette does not need indices.) This is 
one of the simplest models with non-trivial topological content. 

In order to define univocally the action in the lattice simulation, we only need the U ntfl gauge 
fields and do not need to keep track of the phases n ^. To get the topological information we 
could think naively of extracting the $ n phases from the fundamental plaquettes as 

exp(i^) := U n . (2) 

which is ill defined. The ambiguity in the definition of these $ n would give rise to the existence 
of artificial instanton configurations. We will see shortly that if we keep the phases 9 n ^ , a more 
careful definition is possible. 
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In the continuum theory the topological charge in a region v may be computed easily from 
the field strength tensor F^ v via 




(3) 



Note that it is usually understood to take v to be the whole space-time in which case, assuming 
a finite action, one gets an integer winding number. It is not an integer if one only assumes a 
finite action density, so there is no difference from the lattice case with free boundary conditions. 
The charge Q v may be written as an integral over the boundary of v 



with being the normalized tangent vector ofdv. From eq. (4) it is easily seen that for periodic 
boundary conditions in space and time the topological charge of the whole space-time will always 
be zero. 

There are several possibilities in order to define the lattice version of the topological charge 
density and from it Q v . The field theoretical definition for Q v , based on the fact that an expansion 
of the plaquette to lowest order in a gives the field strength[21], is 



and is clearly free of any ambiguities as it involves the plaquette in a well defined manner. 
The disadvantage of this definition is that eq.(4) gets spoiled for finite lattice spacings. Notice 
that this comes from the fact that one does not have strictly a density in the sense that if 
we take double the lattice spacing we do not get exactly the same Q. Therefore, finite action 
configurations would not have exactly integer (or zero, for periodic boundary conditions) winding 
number, but it can get closer for finer lattices. 

We can avoid this problem by defining a geometrical charge density [18, 14], which is the 
simplest example of the fibre bundle construction used in the non-abelian case [19, 10, 20], which 
is in this case just 



This may violate eq.(4) by an integer amount in some cases. In particular you may change 
the winding number by a continuous and local change in the fields. Take for example a periodic 
lattice (Q = 0) with charge on one plaquette close to 1/2, so that if one changes slightly one 
link, the charge on that plaquette will jump to —1/2, the remaining plaquettes staying almost 
the same, thus obtaining Q = — 1 for the whole lattice. These lattice artifacts can be avoided 
for fine enough lattices as we will show. 

In order to keep track of these ambiguities we will directly work with the phases 9 n ^ of the 
links U ntfl as the fundamental variables, from which 9 n is given in obvious notation as the sum 
n = 8n,£: + 8n+£:,T ~ 8n+T,£: ~ n ,r- We then define Q v to be 




(4) 




(5) 




(6) 




(7) 
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where the brackets shift 9 n by an integer multiple of 2ir into the interval ( — vr, it]. This definition 
is clearly equivalent to (6), however, using 9 n as the fundamental variable we are now able to 
locate lattice artifacts (called dislocations), namely places for which \9 n \ > it. Note that in our 
Metropolis algorithm we only update the 9 n ^ by small changes from the old values, otherwise 
we could trivially generate dislocations by adding 2ir to a link, in a plaquette with already low 
action. (These trivial dislocations come as neighboring pairs which do not affect the physics as 
we will see in sec. 5.1 ). The experience has shown us that if we choose parameters corresponding 
to a fine enough lattice, so that the physical instantons are much larger than a, we can avoid 
obtaining any dislocations at all, over the whole Monte Carlo run. It would be interesting to 
know for the non-abelian case, where a characteristic instanton size is only obtained due to 
renormalizaton, if one could also avoid completely dislocations. 

In any case, for non-abelian gauge theories, the geometrical approach to topological charge 
is far more involved, and reasonable criteria to avoid (or exploit) the existence of dislocations 
would be useful, even in order to be able to use safely the simpler field theoretical definition. 
To this end we will consider in our model both definitions for Q v and will show in sec. 5.2 that 
we arrive essentially at the same results. Non-existence of artifacts will be made (extremely) 
plausible for the field theoretical definition and will be controlled by looking at the full 9 n before 
taking the "mod" for the geometrical definition. 

In general an instanton is an extended object and in order to identify this object one has to 
compute the topological charge in a certain region in space-time. This region has to be large 
enough to contain at least one instanton and it should not be larger than half the full lattice. 
Analyzing for example the topological charge of the whole lattice would give a nontrivial result 
only in the presence of dislocations. 

An alternative used by many authors is to identify topological excitations from the local 
topological density. Usually short range fluctuations are much larger than the signal and the 
noise has to be removed by (non-linear) cooling procedures, which affect the excitations too, 
probably giving misleading results. 

One of the physically most relevant quantities to extract is the mean density of instantons, 
p. This can be obtained easily, as we will see, from the charge susceptibility < Q\ > on finite 
volumes v, while for the whole V it is again without dislocations. 

3 Choice of lattice parameters 

In order to be able to extract physical information relevant in the continuum from lattice sim- 
ulations we need that all relevant correlation lenghts be much larger than the lattice cutoff a. 
As for this 2-D model we do not expect a second order phase transition at finite parameters, 
we can at best hope to find a rapid crossover parameter region where the correlations are long 
range, in order to identify extended topological objects without having lattice artifacts. Once 
such a parameter region has been found, we have effectively a fine lattice spacing a and we can 
make sure that all values of 9 n stay small during the whole simulation, without producing any 
dislocations. 

With this aim we analyze the correlation functions of the following three operators [22] on a 
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lattice of size L x X L T with spatial and temporal coordinates x and r respectively. 

Os(T) = Re [J2 %,r)U( x ,r),^( X ,r) + J (8) 

Ov(r) = Ira (£ ^, t) ^, t ) J (9) 

\X,fJ, / 

0"(r) = £ exp^a) Jm^) (10) 

X 

corresponding to the scalar and vector mesons and the photon with momentum kj = 2irj/L x . 
We have to take the connected correlation function in cases where there is a non-zero vacuum 
expectation value < 0{ >. Note that for the 7 we consider higher momentum states which 
couple to the real photon [23], while the O^ =0 operator which corresponds to the field theoretic 
definition for the topological charge in 2-D, gives us almost zero gap energy corresponding to 
the almost degenerate vacuum states which are accessible in the presence of instantons. Even if 
such a picture seems to hold qualitatively ( the correlation of O^ =0 reaches zero mass when the 
coupling k gets low enough so as to start producing instantons), it would be interesting to study 
this more carefully in the future. 

We are not going to give a detailed discussion of the phase structure. Instead we merely 
state that there is a rapid crossover, between regions similar to the Coulomb and Higgs phases in 
higher D. The masses m s , m v and m 1 (k ^ 0) all have a minimum in this region while m 1 (k = 0) 
which is almost degenerate with m v and m 1 (k ^ 0) (right dispersion relation here) at high 
k, drops to zero when we lower k across the region and below. The crossover roughly takes 
place at a set of parameter pairs (A,k) laying on a straight line between the points (.1, .35) and 
(.5, .45). This line of minima is not very sensitive to the gauge coupling (3. As in this paper 
we are interested in a very clear distinction between lattice artifacts and extended (physical) 
instantons, a large value of (3 = 10. has been chosen, corresponding to very fine lattices where 
dislocations practically do not appear anymore. In sec. 5, for the purpose of studying the effect 
of dislocations, we will reduce the gauge coupling to (3 = 6.3. Most simulations have been done 
then around the "scaling region", at (A,k) = (.2, .37), for which we have obtained the masses, 
m s ~ .65, to v ~ .57, m 1 (k ^ 0) ~ .60 and m 1 (k = 0) ~ .02 in lattice units. We are in the region 
where this m 1 (k = 0) drops almost to zero signalling the appearance of instantons as we lower 
k. We will see in the next section that with these parameters we get fairly extended instantons 
with a sizable topological mean density. 



4 Identifying Instantons 

For a clear isolation of instantons it is useful to consider first the case of small instanton den- 
sity, which can be obtained by raising the temperature [24] to high T. Thus we will choose 
narrow lattices in the temporal direction. On the other hand being close to the continuum limit 
implies large instantons in lattice units, and if we want to allow the system to create instanton- 
antiinstanton pairs we need to have large lattices in spatial direction. Consequently we will use 
extremely non-quadratic lattices typically being of the order of 6 X 160. 

In general the observed quantity will be the topological charge Q v in a rectangular volume 
with size L x X L T and with center (x, r). However, in the high temperature regime you would not 
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Figure 1: Inst.-antiinst. pairs are created(destroyed) in Monte Carlo time as seen from the 
sudden rise(fall) in spread in Q v at some two positions, (v = 20 X 4) 




Figure 2: Same run showing the position where the charge Q v is max.(min.) and greater than 
1.0, indicating presence of inst.(antiinst.) there. 
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find strong structural dependence in the temporal direction and we will take L T = N T throughout 
this paper, in which case Q v doesn't depend on r. This has the added nice feature of decreasing 
radically the ultraviolet fluctuations, which grow as the square root of the surface, in this case 
just 2N T ( due to periodicity, the time "boundary" phases cancel). The L x has to be somewhat 
large than an instanton and we take values between 20 and 160 lattice spacings. 

We perfom the simulation with a Metropolis algorithm using a multi-hit updating with 10 hits 
per update with small increments in the phases so as to avoid artificial dislocations. Depending 
on the observed quantity the number of iterations is in the range of 10 4 to 10 6 . We want to 
emphasize that if not otherwise stated the phases 9 n stay safely small (6 n <C vr) for all positions 
and iterations. Thus we are free of lattice artifacts. 

Doing now a Monte Carlo simulation we observe that the system makes rapid jumps from 
a vacuum configuration (including fluctuations) and metastable states consisting of instanton- 
antiinstanton pairs. Depending on the parameters two or more pairs or sometimes even instan- 
tons with charge ±2 are distinguishable. 

To demonstrate this we show in Fig.l for a 4 X 160 lattice the maximal (minimal) measured 
Q v value for any position x, on volumes v (4 X 20). We plot only when \Q\ > 1.0 and averaged 
over the last 10 runs ( just for visualization purposes, only for this plot). We observe in this 
histogram clear jumps for the values of Q v by almost one unit. We interpret this as creation 
or annihilation of an instanton- antiinstanton pair with topological charge ±1. The jump looks 
indeed to be less than one unit because the probability to have a large positive fluctuation on 
top of the instanton charge in a volume v containing it, is much smaller than the probability to 
have the fluctuation somewhere else. 

The dynamics of the pair in Monte-Carlo time is visualized in Fig. 2. We have plotted for 
the same simulation the location of the maximal (minimal) Q v as a function of Monte-Carlo 
time T. For T ~ 6300 a short lived pair appears at location x ~ 80, then there is a period 
with no instantons until at T ~ 8500 a new pair is created at x ~ 120. These excitations move 
fairly independent and randomly through the lattice, even winding around the lattice and finally 
annihilate each other at x ~ 50. The mean size of the (anti)instantons, Li, can be first estimated 
crudely, and will be shown to range between 15 and 30 for most of our simulations. 

For a quantitative analysis of the data we will evaluate the probability distribution P V (Q) to 
get a certain charge Q in a sub-volume v of V, which will be discretized into classes of width 
AQ = .1. For pure fluctuations the distribution is expected to be gaussian around Q = 0. 
As we see from Fig. 3 and Fig. 4 the true distribution is far from being gaussian. Instead we 
observe clear peaks around integer values of Q caused by the (anti)instanton excitations with 
integer topological charge. The peaks are not precisely at integer values due to the overlap with 
gaussian tails caused by the fluctuations and due to the finite size of the instantons, as we will 
show. 

We have fitted the probability distribution of Q with a simple model where we assumed 
(almost) non-interacting instantons and antiinstantons with topological charge ±1 and with a 
uniform charge density on its size Li X N T . 

With p being the average density of instantons (or antiinstanstons), which we want ultimately 
to determine, the probability p n ,m(pV) to have n instantons and m antiinstantons in the whole 
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Top. charge distr. for loops 6*160 in 6*320 latt. 
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Figure 3: Probability distribution to find a top. charge Q in volumes v = 160 X 6, half of V. 
Note the agreement with the data over 3 orders of magnitude. The enhancements at integers 
correspond to multi-instantons fully in v. 




Figure 4: Same as above but with only half the volumes v over which the charge is measured. 
Now the probability is lower to have multi-instantons in v, thus less pronounced peeks than 
above. 
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volume V is equal to the product of two poissonian distributions [15] 



Pn,m(pV) = : ;— exp(-2pV) . 

n\ ml 




(11) 



Taking into account that the boundary conditions enforce equal numbers of instantons and 
antiinstantons the probability p(n) to have n pairs is now given by the diagonal elements of (11) 



where N is the normalization. 

With p n being the probability weight to have n pairs, we have to ask now, in how many 
ways can we get a certain charge Q in the sub-volume v, by displacing the n extended instantons 
and n antiinstantons to all possible locations in V. This number of possibilities to get some Q, 
normalized so as to get with certainty any conceivable one in this sector, gives us the charge 
probability distribution for n pairs, fn(Q). The total charge probability distribution, P V (Q), is 
then just 



For pointlike instantons, this distribution would clearly only give sharp peaks at integer Q's as 
the topological charge would be in or out of v. For extended ones, the charge can be partially in 
v, smearing the peaks around the integers. One can solve analytically f^iQ) for n = 1, (giving 
peaks at and ±1 and a plateau in between), but for higher ra's one has to solve it numerically. 

Up to this point, we have then obtained a charge probability distribution assuming a perfect 
gas of (anti)instantons with a uniform charge density on their size vj. In reality, we have to 
consider also the short range fluctuations which fortunately only grow like ^2N T (no v depen- 
dence at constant temperature). This can be done by assuming that any "exact" charge in v, 
for defining P V (Q), can get smeared by the same gaussian distribution with a certain width a, 
representing the fluctuations. Finally we can include a very primitive, small short range inter- 
action, telling us that instanton- antiinstantons annihilate if they get too close together (< R), 
by just shifting the class in these close encounters to class 0, when computing f^(Q). 

In Fig. 3 and Fig. 4 we have included the best fit of our model to the data (obtained with the 
couplings chosen in last section, corresponding to a fine lattice). We find, using our model for 
a full volume of V = 320 X 6 (solved up to 3 instanton- antiinstanton pairs), a good fit for both 
sub-volumes v, with the parameters in lattice units p ~ .0013, a ~ .27, Lj ~ 28 and R ~ 4. The 
fact of having used only 3 pairs to calculate the fit, makes an overestimation for p which should 
shift this value down to p ~ .0011 in the exact case, close to the value that we will get from the 
susceptibility. 

In the following we will calculate the second moment of P V (Q), < Q\ >, which will be easier 
to calculate than the full P and will give us directly p, for the case of non-interacting instantons. 

The contribution p I,A {Q) of a single instanton (antiinstanton) in a volume V with uniform 
distributed position to the total probability distribution of Q measured in v (v > Lj) includes 
the probabilities to be completely in v (p; n = v ~ I y Nt ), completely out of v (p Q ut = v ~ L \^ t ~ v ) 
and to lie on the surface (p sur f = 2L y Nt ) , finally resulting in 



Pn(pV) 



(pV) 2n 1 
nl 2 N ' 






n 




(14) 
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To this distribution of excitations we have added independent gaussian quantum fluctuations, 
denoted by p G (Q) of the topological charge with width a, which is in general a function of v. 
The expectation value of Q 2 can be calculated exactly to give (using p\Q) = p A { — Q)) 

oo poo rl 2n n 2n 

<Ql> = Y,Pn(pv) d Q I[dQ j (Q + Y l Qj-Qn +j ) 2 P G (Q)'[[p I (Q j ) 

ri J — OO JO - ^ 

= * 2 + 2(£(l-£)-^)£»P»(^)- (15) 

n 

The sum in the last expression may be approximated to J2 n n Pn(pV) ~ pV for pV 1. 

We did several Monte Carlo runs for volumes at various temperatures T, 
V = 320 X (4,6,8,10,12) , and for various subvolumes v = (60,80,120,160) X N T in order to 
check eq.(15). We found indeed, for this range of v for a fixed N T = 1/T, that we could extract 
a fairly stable p(N T ). There is a temperature dependence, as expected [24] from the squeezing 
of instantons at high T, obtaining respectively: 

/o(4, 6, 8, 10, 12)= (0.4(2), 1.1(2), 1.4(2), 1.4(1), 1.5(2)) X 10" 3 . This instanton density p grows 
fast as we decrease the temperature and starts stabilizing at low temperatures as wanted. Note 
that the density for N T = 6 obtained here agrees well with the one extracted from P V (Q). 

Since topological charge in a volume v may be expressed as an integral over the boundary of v, 
statistically independent fluctuations within v do not affect Q 2 and hence a 2 will be proportional 
to the surface of v for sufficiently large v. In our case dv consists of two rings with a fixed length 
running around the torus. If the separation v x of the two lines is sufficiently large the fluctuation 
parts of the corresponding loop integrals, denoted by < K(x) > and < K(x + v x ) >, are 
independent and a does not depend on v x . This formal result is confirmed by our numerical 
calculations. 

In contrast, the instanton contribution to < Q 2 > has (apart from finite size effects, which 
is only effective on the surface) a different statistical behavior because the instanton parts of 
the loop integrals can never be independent. (If they were < Q 2 > and thus the action density 
would be infinite. In fact, formally speaking, the loop integrals are maximally correlated for any 
v x since we have for the covariance 

< K{x)K{x + v x ) > _ _ 1 < Q 2 {v) > _ 

V< K(x) 2 >< K(x + v x ) 2 > 2 ^< K(x) 2 >< K(x + v x ) 2 > 

which is merely a consequence of the fact that the number of vacua is infinite.) 

The different statistical behavior of the instanton and the fluctuation contribution can be 
used to distinguish the two. We believe that this is in general true as long as you can express 
topological charge in a volume v in terms of the boundary of v. 

5 Beyond optimal conditions 

5.1 The role of dislocations 

So far the simulations were done for parameters for which in practice any configuration was 
absolutely free of dislocations. Dislocations can occur when after an update a plaquette angle 
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n 





1 


2 


3 


with dislocations 


1. 


.51±.04 


.11±.015 


.005±.002 


without dislocations 


1. 


.53±.02 


.10±.01 


.007±.003 



Table 1: Relative weights of sectors with charge n 

9 n starts fluctuating around 2irn with a dislocation charge n ^ 0. In general due to the modulus 
in eq.(7) for N dislocations with charges n 4 - the topological charge of the whole lattice is 

N 

Q tot = -]>>«. ( 16 ) 

8 = 1 

A proper choice of parameters alone is not sufficient to avoid dislocations. A hot start of the 
Monte-Carlo simulation with random link variables will always produce dislocations being more 
and more stable the closer one comes to the continuum limit. 

Since you can move a dislocation by adding an integer multiple of 2ir to an angle of a 
link variable, which affects neither the action nor the topological density, it is clear that the 
system does not 'know' about the position of a dislocation. Furthermore, since dislocations of 
opposite sign (and equal magnitude) can be moved through the lattice to cancel each other, the 
only relevant quantity can be the net dislocation charge, which merely forces the 'real' total 
topological charge to assume the value Q tot and thus effectively acts like a different boundary 
condition (costing some action). When doing a Monte-Carlo simulation the configurations may 
now be divided into a discrete set of sectors with fixed Q tot . Now we will address the question 
of whether a lattice of volume v with inclusion of dislocations, behaves like a subsystem (with 
the same volume) of a larger lattice V without dislocations. 

When doing a careful analysis we must first of all make sure that we are in a parameter 
range where the jumps to different sectors are rare, for otherwise these jumps itself would spoil 
thermal equilibrium of each individual sector. We achieved this by choosing (3 = 6.3 without 
changing A and k. For this (3 a jump occurs approximately every 300 iterations. 

To check numerically whether only the net dislocation charge can affect the system we partly 
modified the lagrangian to prohibit occurences of dislocations without changing the continuum 
limit. To be specific, we added a term c • 0(|# n | — 7r), with c being a large constant, to the 
lagrangian density. We now compared results with the modified lagrangian, to the Q tot = 
sector of the original lagrangian, choosing all parameters identical and we found that all observed 
quantities including the probability distribution of Q come out identically. 

Furthermore, though the absolute amount of dislocation charges always increases (like a 
random walk), the distribution of the net dislocation charge, denoted by P(n), converges with 
Monte-Carlo time and is (relative to the former) highly peaked around Q tot = 0. We compared 
the distribution of Q tot , obtained in a 40 X 4 lattice, with results from the modified version where 
we have measured the distribution P V (Q) in a 40 X 4 subvolume of a 320 X 4 lattice. For the 
latter the separation into sectors of charge n was simply done by integrating P V (Q) over intervals 
[n-\,n+ \}. 

The resulting probability distribution of the sectors is shown in table 1. To make the relative 
weights of the sectors obvious we have normalized both distributions to P(0) = 1. The agreement 
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Figure 5: Both definitions give similar curves for fine lattices, the field theoretical one losing 
some information at large Q. 

of the two methods, within the errors in Table 1, confirm our hypothesis that dislocations simulate 
a larger volume in which the lattice can be thought to be embedded. This especially means that 
the dislocations seem to come with the right entropy. Therefore, the more usual way to measure 
the topological charge, based on using the dislocations to get non-zero charge seems fine, at 
least when the physical instantons are much larger than the lattice spacing. Clearly when the 
instantons are of the size of a, their action gets distorted and small, as they are not distinguishable 
from the dislocations anymore. Then the picture above loses meaning, but we were hopelessly 
out of the continuum scaling region. 

5.2 The field theoretical definition of topological charge 

For non-abelian theories it is difficult to implement the geometrical definition of topological 
charge [19, 20, 10]. A far more simple definition is given by the field theoretical one [21] in which 
the topological density (in general tr(FF*)) is expressed in terms of the plaquettes with the 
correct continuum limit. 

The most severe problem with the field theoretical definition in eq.(5) is the (in contrast to 
dislocations uncontrollable) violation of eq.(4). Fluctuations within a volume v do not cancel 
exactly, as happens with the geometrical definition in eq.(6). As a consequence, the fluctuation 
part of Q measured in v will, like the instantons themselves, have a contribution proportional to 
y/v. The different statistical behavior of instanton and fluctuation contributions to Q can not 
be used to distinguish the two and the deviation from the gaussian distribution of Q, the peaks 
at integer values, will be washed out, the larger one chooses v. 

However, the difference of the two definitions of topological charge J2n( s ' m (@n) — [&n]) ' 1S only 
of order [0 n ] 3 and there is hope that for small fluctuations the peaks are still visible. In our 
case this is in fact true for moderate v as is shown in Fig. 5. Here we see good agreement of the 
probability distributions. In the non-abelian case though, one expects a poorer agreement [21], 
as the deviations from the geometrical case do not decrease for finer lattices. 
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6 Conclusion 



We have seen that, at least in the two dimensional case, one can obtain extended physical 
instantons in a periodic box without any lattice artifacts. These instantons can be identified 
and their mean density p and size measured, by considering the charge Q v in sub-volumes v of 
the full lattice. The probability distribution to get some topological charge in v, P V (Q), can 
be obtained by lattice simulations and gives clear enhancements at integer Q. A simple model 
of almost non-interacting extended (anti)instantons fits nicely the Monte Carlo results. From a 
volume analysis of the susceptibility < Q\ > we can easily extract the mean instanton density 
p and find that it grows as we decrease the temperature until reaching some value at zero T, as 
expected. 

Furthermore, we have seen that if one changes the couplings so as to allow some dislocations, 
still for fine enough lattices, the system acts as if we had changing boundary conditions for some 
number of physical instantons, having practically the same probability distribution as in the case 
above for instantons in a sub-volume v of the same size and very large V. This is the more usual 
way to get topological charge distributions, but it was interesting to check that it really gives 
similar answers to the proper distribution without lattice artifacts. 

This whole idea of measuring the topological charge on subvolumes of the periodic lattice, 
with v much larger than the ultraviolet fluctuations correlation lenghts and even larger than the 
extended instantons and smaller than V/2 , could be similarly tried out for non-abelian theories 
in higher dimensions. Clearly, the geometrical definition for topological charge is preferable at 
least as control. One would need though, fairly large average instantons sizes, in order to avoid 
dislocations. The larger dispersion in the sizes of instantons, compared to our 2-D case, might 
be the biggest problem. Improved actions certainly would help. The lattice would only need to 
be elongated in one direction, in order to decrease the fluctuations. Finally, the field theoretical 
definition could be useful at intermediate volumes for fine enough lattices. 
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